Critical phenomena of thick branes in warped spacetimes 



Institute of Cosmology and Gravitation, 



Antonio Campos 
University of Portsmouth, Portsmouth POl 2EG, United Kingdom 
(Dated: February 1, 2008) 



We have investigated the effects of a generic bulk first-order phase transition on thick Minkowski 
branes in warped geometries. As occurs in Euclidean space, when the system is brought near 
the phase transition an interface separating two ordered phases splits into two interfaces with a 
disordered phase in between. A remarkable and distinctive feature is that the critical temperature 
of the phase transition is lowered due to pure geometrical effects. We have studied a variety of 
critical exponents and the evolution of the transverse-traceless sector of the metric fluctuations. 
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The physics of branes in higher dimensional theories 
has been the subject of renewed interest in the past years, 
mostly, because they provide a novel approach for resolv- 
ing the cosmological constant and the hierarchy problems 
jj], 0. An intriguing feature of this type of scenario is 
the possibility of localizing a massless graviton on the 
brane, and then reproducing effectively four-dimensional 
gravity at large distances, due to the warp geometry of 
the spacetime ||. Furthermore, these models have also 
inspired the idea in which our universe is the result of 
the continuous collisions of branes [|J. A fundamental 
ingredient of this cosmological model is the nucleation 
and splitting of branes. Then, it is interesting to con- 
sider the universal aspects of thick branes p| splitting in 
a warped bulk. In this work we present a field theory de- 
scription for a feasible realization of such a mechanism. 
Using superpotential techniques ^ n], some exact solu- 
tions of the equations of motion for the coupled scalar 
field plus gravity system have been studied in 0. Ana- 
lytical solutions for giving potentials have been described 
in (pi. The probability of tunneling has been considered 
in [p| and some geometrical aspects of brane collisions 
have been pointed out in [|l0| O. 

In this letter, we provide a qualitative description of 
the splitting of thick Minkowski branes due to a first- 
order phase transition in a warped geometry. For this 
purpose we have considered a generic potential for a com- 
plex scalar field that captures the universal aspects of the 
phase transition in the bulk. As the system is brought 
near the critical temperature an interface interpolating 
two bulk phases breaks into two separated interfaces, and 
a layer of a new phase appears between them. This effect 
is known in condensed matter physics as complete wet- 
ting |l3| . It also appears at the deconfinement phase 
transition of SU(3) Yang-Mills theory |l4], [f5| and super- 
symmetric QCD |l6| . Our analysis shows the existence of 
a new effective critical temperature for complete wetting 
due to the effective negative bulk cosmological constant. 
This phenomenon reveals a remarkable relationship be- 
tween geometry and the critical behavior of interfaces. 
We have been able to characterize the universal aspects 
of this critical phenomenon computing the critical expo- 



nents for some relevant quantities. And finally, we have 
also studied the behavior of the metric fluctuations. In 
particular, we have found that the zero mode delocalizes, 
high energy modes are scarcely affected, and for low en- 
ergy modes the phase changes and the amplitude starts 
to grow at the moment of the splitting. It is interesting 
to point out that we are going to use a numerical ap- 
proach and thereby our solutions do not depend on any 
specific superpotential. 

We start with a metric of the form 



ds 2 



e 2Air) r Hj dx i dx j 



dr z 



(1) 



where i = 0,1,2,3, r)ij = (1, — 1, — I, — I), and the warp 
factor A(r) depends only on the coordinate of the extra 
dimension r. This choice of metric preserves Poincare 
invariance in the four-dimensional submanifolds of con- 
stant r. To build up the thick brane we consider an aux- 
iliary complex scalar field $ coupled to gravity through 
the five-dimensional action 

S = J d 4 x dry/\f\(~R+^ v d^*d^ -V "($)- A), 

where R is the scalar curvature and A is a bulk parameter 
that fixes the minimum of the scalar potential V($). In 
flat backgrounds A is a free parameter that can always 
be used to set the minimun of the potential to zero. In 
contrast, this is not always possible if one is concerned 
to find out soliton solutions of the equations of motion in 
curved spacetimes. Since we are particularly interested 
in the universal aspects of first-order phase transitions in 
the bulk we are going to consider the generic scalar field 
potential 



am 



3#) + c|#| 
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where 4>r and <fii are, respectively, the real and the imag- 
inary part of the scalar field <f>. This potential has quar- 
tic self-interactions and respects Z3 symmetry. The ac- 
tion for a scalar field with a potential of this type be- 
longs to the same universality class as the effective ac- 
tion for theories describing the Zjv center symmetry of 
a SU(N) gauge field ]l5| |. So, they can be easily re- 
alized if the model supports a bulk field with a gauge 
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symmetry. To ensure stability the parameter c must 
be positive. For a > the potential has a local min- 
imum at $W = corresponding to a disordered bulk 
phase. For b > and a < 9a c /8 there are three ad- 
ditional degenerate global minima corresponding to or- 
dered bulk phases, say $ (2) = <j> , $ (3) = (-± + i^)(f> , 
and <£> (4) = {-\ - i^)<t>o, where 

3b ( I 8a~\ 
*° = 8c (W^ 9*)' (3) 

and a c = b 2 /4c. These three minima coexist with the 
disordered phase when a = a c . Then fixing b and c 
lead us to interpret a as a bulk temperature parameter 
where a c is the critical temperature. The actual vari- 
ation of a will strongly depend on the particular real- 
ization of the effective theory. But, since we are study- 
ing the universal aspects of the phase transition, we can 
treat this parameter as an independent variable. Note 
that the effective five-dimensional cosmological constant 
is not A but A e // = A + V(a), where V(a) indicates the 
value of the scalar potential at the global minima, i.e. 
V(<&( n )) = cj) 2 (a — b(f) a + ccj> 2 ) with n taking values 2, 3 or 
4. 

The classical equations of motion for planar field con- 
figurations read as follows || , 

4>" R = -4A'<p' R +^-, (4) 

0(f>R 

= -4A% + (5) 
d(pi 

A" = ~(<P'l + <!>?), (6) 
3A' 2 = lfrZ + <tf)-V(*)-A, (7) 

where the prime denotes differentiation with respect to 
the coordinate r. The soliton equations for the flat case 
are recovered by taking the warp factor A(r) to be con- 
stant |l5| (see [|l6| for a supersymmetric generalization). 
Note also that Eqs. (g) and (§) imply that A < \V(a)\ 
and then the effective five-dimensional cosmological con- 
stant A e fj is always negative. We have solved numer- 
ically the system of ordinary differential equations (|])- 
(Q) for (f>R(r), <f>i(r), and A{r) with boundary conditions 
$(-oo) = $< 4 \ $(-oo) = $( 3 ), A(0) = A'(0) = (see 
Fig. |l|). Then, Eq. (Q) has been used in order to fix 
the parameter A for different values of the parameter a 
as the interface configuration is brought near the phase 
transition. Because of the Z3 symmetry there will be 
equivalent solutions for other choices of ordered phases. 

Our results reveals that, equivalently to what hap- 
pens when the background spacetime is Euclidean, the 
ordered-ordered domain wall splits into two separated 
ordered-disordered interfaces, and a layer of disordered 
phase appears and completely wets the two ordered 




-20 -10 10 20 



r 

FIG. 1: Interface profiles for thick domain walls in a warped 
geometry with a Z3-symmetric bulk first-order phase transi- 
tion for b — 2 and c = 1. (a) Away from the phase transition 
a domain wall interpolates two distinct ordered phases, (b) 
Close to the critical temperature the wall splits into two inter- 
faces with an intervining film of disordered phase in between. 

phases as the the phase transition is approached. Al- 
though in both cases this behavior looks quite similar 
there is an important distinctive feature. For a geometry 
with a warp factor complete wetting does not occur for 
the critical temperature a c but for a smaller effective tem- 
perature a*. As we will see, this effect is closely related 
to the fact that asymptotically the spacetime is Anti- 
deSitter. In order to understand why this effective critical 
temperature appears it is necessary to know the behavior 
of the Euclidean energy E(r) = + tfif) — V($) as 

the temperature parameter a increases. It is very easy to 
see, using the equations of motion (|J) and (||), that the 
derivative of this energy along the orthogonal direction 
to the domain wall is E'(r) = -4 A' ((/>'£ +$f). Then, it 
has a minimun at r = 0, -E'(O) = A, and it reaches its 
maximun value at infinity _E(±oo) = | V(a) | . In fact, it 
has the same shape as 3A' 2 but is displaced vertically. On 
the other hand, one also has to take into account that the 
stability of a configuration with three coexisting phases 
implies that an interface of two ordered phases can only 
be completely wet by a disordered phase. This means 
that in order to have a disordered phase between our ini- 
tial ordered phases $(0) — > $W and $'(0) — » when a 
increases. Moreover, E(0) should tend to zero as well as 
A for complete wetting to occur. If the background is 
Euclidean, the energy is a positive constant, E = \V(a)\, 
for each value of the temperature parameter a. So, since 
E is zero for V(a = a c ) — 0, complete wetting happens 
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with a critical temperature a c . However, in the warped 
case and because A is always smaller than | V(a) | , the 
coexistence of the three phases occurs for a critical value 
a* < a c . This inmediately translates into the fact that 
the effective bulk cosmological constant asymptotically 
reaches a constant negative value 



A 



eff 



V(o)|«. < 0. 



(8) 



Then, one could revert the argument and say that the 
existence of a nonzero negative cosmological constant is 
responsable for lowering the critical temperature of com- 
plete wetting. One would expect intuitively that the 
value of a* is very close to the critical value a c . In this 
case, it can be easily seen that A*^j, a Cl and a* are re- 
lated by, 



A 



eff 



a c , 
— \a c 
c 



a*)- 



(9) 



It is worth noticing that the existence of an effective crit- 
ical temperature prevents the global geometry for be- 
coming flat as the phase transition is reached. If the 
contrary would have occurred, it would have been a puz- 
zling feature of brane models. In the example plotted 
in Fig. [I] the critical temperature is a c = 1. In this 
case, we have obtained the effective critical temperature 
a* = 0.9223822740 for A < 10" 7 and an effective bulk 



cosmological constant A 



eff 



-0.08. 



The universal aspects of this critical behavior are char- 
acterized by several critical exponents. In Fig. ^ we have 
plotted the numerical results for the determination of the 
critical exponents of several parameters of the phase tran- 
sition. As one would expect from the discussion above, 
the parameter A is proportional to the deviation from the 
effective critical temperature, (a* — a)/a* oc A. On the 
other hand, the width of the disordered phase layer, r„, 
wetting the two ordered phases diverges logarithmically 
r Q oc — log(A). In addition, for the order parameters at 
the center of the wetting layer we obtain A"(0) oc A, 
4>r(0) oc A 1/2 , and 4>'j(0) oc A 1/2 . In condensed mat- 
ter physics, this type of critical behaviour corresponds to 
systems with short-range interactions |lj| . 

Now, let us consider the metric fluctuations with re- 
spect to this background solution. In the axial gauge 
ha = = /i44 the perturbed metric can be written as 



ds 2 



^ A(r Hmj + h^dx'dx 3 - dr 2 



(10) 



The transverse and traceless (TT) sector of the perturba- 
tions, hij, decouples from the scalar field perturbations 
and satisfies the equation p], 



(d 2 + AA'd r - e- 2A U) hij = 



0. 



(11) 



This equation can be transformed into a Schrodinger-like 
equation by changing to the coordinate dz = e 
introducing the new variable Hij = e 



-ikx e (3A/2] 



A dr and 

hij, 



d 2 m 



k ) m 



(12) 
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FIG. 2: Determination of the critical exponents for a variety 
of order parameters of the phase transition. 



Here, the potential is 



U(z) 



9 A2 



-A 



(2A" + bA'' 



(13) 



where the dot now stands for the derivative with respect 
to the new coordinate z. The zero mode k = has the 
analytical solution H^ — Nije^ A / 2 \ with Nij a normal- 
ization factor. As shown in Fig. away from the phase 
transition the potential in the Schrodinger equation for 
the TT sector of the metric perturbations has the shape 
of a volcano potential. However, after the wall splits off, 
the shape becomes that of a double-well-type potential 
with a depth decreasing as the two walls separate. One 
can observe that the critical behavior delocalizes the zero 
mode, but it does not introduce any gap between the 
zero mode and the continuum spectrum. Furthermore, 
the zero mode is the lowest energy eigenstate and then 
there are no instabilities coming from tachyonic modes. 

In Fig. ||, we have also plotted the wave functions of 
two different types of nonzero modes. In Fig. ||(b), we 
have considered the case of modes with energy larger than 
the maximun of the volcano potential, and in Fig. |](a), 
those with a smaller energy. The first thing to note is 
that the wave functions get broader as the mode has 
lower energy. Modes with larger energy barely get af- 
fected by the splitting of the interfaces. On the contrary, 
the phases of the waves for modes with a small energy 
change at the moment of the nucleation, increasing their 
amplitude afterwards. 

In conclusion, we have investigated the critical behav- 
ior of coexisting interfaces in an asymptotically anti-de 
Sitter bulk. Because a film of disordered phase is formed 
between two ordered phases, this mechanism of complete 
wetting gives a generic field theory description of the 
splitting of a thick domain wall, which is a fundamental 
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ingredient for the brane cosmological models discussed 
in P]. Moreover, we have found that naturally an ef- 
fective critical temperature shows up due purely to the 
global geometry. This is a remarkable result because it 
shows a direct link between geometry and the universal 
aspects of a critical process. Nevertheless, there still are 
a few interesting questions that deserve a deeper investi- 
gation. For instance, what is the behavior of the non-TT 
and the scalar field perturbation sectors during the split- 
ting of the thick brane? How does the critical behavior 
found here change if we consider de Sitter or anti-de Sit- 
ter branes instead of Poincare branes? Is, then, multilo- 
calization of the zero mode possible? Which is the effect 
of a bulk scalar field? 

The author is deeply in debt to C. F. Sopuerta and 
very thankful to his former collaborators K. Holland and 
U.-J. Wiese. This work has been supported by the Marie 
Curie Fellowship HPMF-CT- 1999-00 158. 



FIG. 3: Shapes for the volcano potential (a) and the wave 
function of the zeromode (b) for the transverse-traceless part 
of small gravitational perturbations. The critical temperature 
is approached from bottom to top in both plots and we have 
taken an arbitrary normalization for the wave function. Note 
that, away from the transition, the zeromode is bound to the 
walls in the fifth extra dimension but, as the walls peel off 
and pull apart, it is delocalized. 
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FIG. 4: Wave functions for nonzero modes of metric fluctu- 
ations with an arbitrary normalization, (a) Modes with en- 
ergy smaller than the maximun of the volcano potential, and 
(b) modes with larger energy. The plots represent different 
stages of the phase transition in the following order: solid line, 
dashed line, dot line, and dashed-dot line. Being the latest 
the closest to the critical temperature. 
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